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Elliptic flow and incomplete equilibration at RHIC 
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We argue that RHIC data, in particular those on the anisotropic flow coefficients V2 and V4, suggest 
that the matter produced in the early stages of nucleus-nucleus collisions is incompletely thermalized. 
We interpret the parameter (1/ ' S)(dN/dy), where S is the transverse area of the collision zone and 
dN/dy the multiplicity density, as an indicator of the number of collisions per particle at the time 
when elliptic flow is established, and hence as a measure of the degree of equilibration. This number 
serves as a control parameter which can be varied experimentally by changing the system size, the 
centrality or the beam energy. We provide predictions for Cu-Cu collisions at RHIC as well as for 
Pb-Pb collisions at the LHC. 
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The observation of azimuthal anisotropy in the produc- 
tion of particles in ultra-relativistic nucleus- nucleus col- 
lisions, especially the so-called elliptic flow vi Q, is one 
of the highlights of the RHIC heavy ion program 0, Q . 
The phenomenon, first identified in this regime at the 
CERN SPS [5|, reflects the anisotropy of the region of 
overlap of the nuclei and is a direct consequence of the 
reinteractions between the produced particles. In the 
limit where the collisions are frequent enough to drive 
the system quickly to local equilibrium, fluid dynamics 
provides an intuitive physical explanation for the origin 
of anisotropic flow: particles tend to go in the direction 
of the strongest pressure gradients, hence preferably in 
the collision plane [f| . Of course, local equilibrium is not 
a necessary condition for elliptic flow, but it is commonly 
accepted that deviation from equilibrium can only reduce 
the magnitude of the effect. 

In fact, the main characteristics of the observed elliptic 
flow are reasonably well described by ideal fluid dynam- 
ics 0, while requiring unreasonably large cross sections 
in transport models The ability of the former to re- 

produce both the elliptic flow and single-particle spectra 
for measured hadrons with px 2 GeV/c near midrapid- 
ity in minimum-bias collisions is considered a significant 
finding at RHIC; by contrast, at the SPS a simultane- 
ous fit of both observables appears impossible. The de- 
pendence of the flow pattern on hadron masses further 
supports the hydrodynamical picture. Very strong con- 
clusions have been drawn from this apparent success : 
one argues that local equilibrium has to be established on 
short time scales (~ 0.6 fm/c 6]) while viscosity should 
be negligible, suggesting a "perfect fluid" behaviour for 
the created matter [10). 

In this Letter we would like to question these conclu- 
sions. On the one hand, the short time scale for equi- 
libration is difficult to account for microscopically 11] 
(although it has been argued recently that plasma in- 
stabilities may provide a mechanism for fast thermaliza- 



tion; see, e.g., Refs. 0,0). On the other hand, several 
features of the data clearly signal the breakdown of the 
hydrodynamical description and are more naturally un- 
derstood if the constraint of local equilibrium is relaxed. 
Thus, the point of view that we shall adopt here is that 
matter produced at RHIC is not fully equilibrated, and 
we shall explore the consequences of such an assumption. 
As we shall see, this leads to simple predictions that can 
be easily tested. 

This Letter is organized as follows: We first recall the 
essential features of the elliptic flow using an ideal hydro- 
dynamic picture. We then show that incomplete thermal- 
ization leads to specific deviations from hydrodynamical 
behaviour: these concern in particular the dependence 
of moments «2 and of the azimuthal distribution on 
the system size and the collision centrality. Note that 
throughout this paper, we shall only consider the average 
values of i>2 and over all particles at a given rapidity; 
the effects of partial thermalization on the elliptic flow of 
identified particles, in particular its p-r-dependence, have 
been discussed in Ref. pL4| . 

Elliptic flow originates from the anisotropy of the ini- 
tial matter distribution. In hydrodynamics, the depen- 
dences of the elliptic flow on the system size and the 
centrality are essentially determined by the spatial ec- 
centricity, defined as 



w 2 



(y 2 +. 



(1) 



where x and y are coordinates in the plane perpendicu- 
lar to the collision axis (with the ^-direction in the col- 
lision plane, and the origin midway between the centers 
of the two nuclei). Angular brackets (•) denote an aver- 
age weighted with the initial entropy density s{x,y) (at 
a given rapidity). 

Elliptic flow develops gradually in the system as it 
evolves. If the speed of sound c s is constant, the natural 
time scale is of the order of R/c s , where R is a measure 
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FIG. 1: (Color online) i>2/e as a function of c s (t — to)/R 
for various impact parameter values of Au-Au collisions at 
y/s^N = 200 GeV with c s = c/sfl. 



of the transverse size of the system: the anisotropy of the 
momentum distributions can be fully achieved only once 
all parts of the system are "informed" about the initial 
spatial anisotropy, and that takes a tim e of the order of 
R/c s . We define R through l/R = y/l/{x 2 ) + l/(y 2 ) 
(since flow is an effect of pressure gradients, R is a more 
natural choice than, e.g., the rms radius). 

Hydrodynamical results are presented in Fig. ^ which 
displays the time evolution of elliptic flow (i.e., the value 
of v 2 that one would observe if the system was to de- 
couple at time t) for various centralities. The proce- 
dure is detailed in Ref. 5]. There, v 2 is defined as 
v 2 = (Px~Py)/(Px+Py)i where average values are taken 
over all particles. This yields values of v 2 that are typi- 
cally a factor 2 larger than those obtained with the more 
conventional definition, v 2 = ((p x ~p1)/{p x +Py))- The 
longitudinal expansion is assumed boost-invariant; the 
hydrodynamical evolution starts at a time to = 0.6 fm/c 
after the collision. We checked that the results are inde- 
pendent of the value of to as long as to is much smaller 
than R/c s . In line with the discussion above, we have 
plotted v 2 divided by e, and the elapsed time t — to di- 
vided by the characteristic time R/c s . This results in an 
almost perfect scaling for a large range of impact param- 
eters (from b — 2 to b = 12 fm) for which the eccentricity 
varies by more than one order of magnitude and R by 
a factor 2 (see Table [I] below) . Note, in particular, that 
the final value of v 2 is independent of the system size (R) 
for a given shape (e). This is a consequence of the scale 
invariance of ideal fluid dynamics. 

The magnitude of v 2 also depends on the fluid proper- 
ties, in particular on the speed of sound c s . In order to 
illustrate this dependence, we have repeated the calcula- 
tion for various values of c s . The results are presented 
in Fig. |2 When c s is large enough (c s ^ 0.3), v 2 is 
proportional to c s . This no longer holds when c s be- 
comes small (c s ^ 0.2): in this "nonrelativistic" regime, 
the expansion of the system is entirely controlled by the 
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FIG. 2: (Color online) V2/e as a function of c s (t — to)/R for 
various values of c s in the case of Au-Au collisions at ^/sjvjv = 
200 GeV with impact parameter 6 = 8 fm. 



dimensionless parameter c s (t — to)/ R: we have checked 
that v 2 /e becomes a universal curve for c s < 0.1. Note 
that the results plotted in Fig. |21 have been obtained by 
assuming that c s is constant throughout the evolution. 
In the real world, c s varies. In particular, if the system 
enters a regime where the equation of state is soft, i.e., 
where the speed of sound is small, the resulting v 2 may 
be significantly reduced |l5|. 

In summary, we have seen that within ideal fluid dy- 
namics the final value of the elliptic flow v 2 is propor- 
tional to the initial spatial eccentricity e 0, is indepen- 
dent of the system size R, grows with the speed of sound 
c s , and has a characteristic build-up time ~ R/c s . Note 
that this time scale for the build-up of the elliptic flow 
is of the same order of magnitude (although somewhat 
larger) in transport calculations (see, e.g., Refs. |8Hlq). 

We start now exploring the consequences of incomplete 
equilibration. To do so, we characterize the degree of 
thermalization by a dimensionless parameter, the Knud- 
sen number K By definition, K^ 1 is the typical 

number of collisions per particle. Local thermal equilib- 
rium is achieved if K _1 ^> 1. We are going to show that 
K^ 1 can be determined from the data, as it is propor- 
tional to (1 / S) dN I dy , where dN/dy is the total multi- 
plicity density, while S = 1~K\J (x 2 ) (y 2 ) is a measure of 
the transverse area of the collision zone (with this defini- 
tion of S, which is larger by a factor 2 than that adopted 
in Ref. {1/ S) dN/dy is the maximum value of the 
density for a Gaussian density distribution). 

The relevant length and time scales for elliptic flow are 
R and R/c s , respectively. Therefore, the typical num- 
ber of collisions per particle is given by R/X, where A is 
the mean free path at time R/c s . This mean free path 
A depends on the particle density n: 1/A = an, with 
a a cross section characterizing the interactions among 
the produced particles. The particle density, in turn, de- 
pends on time. We assume that the total particle number 
is conserved throughout the evolution: this is justified 
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TABLE I: Parameters for Au-Au collisions at y / ijvjv 
200 GeV. See text for details. 



b (fm) 


e 


R (fm) 


dN 
dy 


i *£ (rnb" 1 ) 


n {2,\ fm -3 








2.07 


1050 


1.95 


5.4 


2 


0.033 


2.02 


975 


1.90 


5.4 


4 


0.115 


1.89 


790 


1.77 


5.5 


6 


0.215 


1.68 


562 


1.55 


5.3 


8 


0.315 


1.45 


344 


1.23 


4.9 


10 


0.398 


1.22 


167 


0.82 


3.8 


12 


0.433 


1.04 


55 


0.37 


2.0 



TABLE II: Parameters for Cu-Cu collisions at ^/sjvjv 
200 GeV 



6 (fm) 


e 


R (fm) 


dN 
dy 




n(£) (fm" 3 ) 








1.42 


275 


1.09 


4.5 


2 


0.043 


1.36 


240 


1.02 


4.3 


4 


0.141 


1.23 


159 


0.83 


3.9 


5.5 


0.216 


1.10 


95 


0.61 


3.2 


6 


0.237 


1.06 


77 


0.53 


2.9 


8 


0.265 


0.93 


22 


0.20 


1.2 



by the observation that n is proportional to the entropy 
density, and that entropy is conserved. This allows us to 
estimate the particle number density at time r from the 
relation 



ctti(t ) 



1 dN 
S dy : 



(2) 



valid for times r R/c s , i.e., as long as the transverse 
size of the system does not vary significantly. One even- 
tually obtains 

K = - = Ran — ) = -x^ ■ (3) 

A \ c s J b dy c 

Tables [I] and |H] provide numerical values of 
(1/5) dN/dy for two colliding systems, Au-Au and Cu- 
Cu, at the same center-of-mass energy ^snn — 200 GeV. 
They are given in mb -1 , making the conversion into a 
Knudsen number easy once the cross section is known 
(up to a factor c s /c). With a typical partonic cross 
section of 3 mb, and a speed of sound c s ~ c/%/3, one 
thus finds K^ 1 ~ 2 for a semi-central Au-Au collision. 
This is not much larger than unity, and the ideal fluid 
limit probably requires larger cross sections, as already 
inferred from transport calculations. The Knudsen num- 
ber can also be estimated from the viscosity [l8j |. with 
similar conclusions. 

Incomplete equilibration breaks the scale invariance of 
ideal fluid dynamics: V2 depends on the system size R 
through K^ 1 . As K , increases, the magnitude of V2 
grows linearly — the larger the number of collisions, the 
larger the momentum anisotropy — , eventually saturat- 
ing when the system reaches local equilibrium. To fix 
ideas, one may have in mind the following simple for- 
mula, which exhibits the correct qualitative behaviour: 



t'2 



hydro 



K- 



K~ 



(4) 



where D2 ydro is the value of «2 obtained from hydrody- 
namics, i.e., corresponding to the large K^ 1 limit, and 
Kq 1 is a number of order unity, whose precise value can 
only be determined through an explicit transport calcula- 
tion. Referring to the transport calculations in Ref. ^(| , 
one finds only a modest increase of v 2 (by 40%) when 
a increases from 3 to 10 mb, K^ 1 from 2 to 7, which 
corresponds to K^ 1 ~ 1.5. 

We are now in position to examine whether the matter 
produced at RHIC is in local equilibrium. To do so, we 
need to vary the value of K~ l . This can be done by 
changing the system size R, while keeping the mean free 
path A constant [see Eq. JSJ]. Since A depends on the 
density (explicitly, and through a potential dependence 
of the cross section on the density), this requires keeping 
the density constant. Now, the last column of Table |U 
shows that the particle density remains approximately 
constant in Au-Au collisions for a large range of impact 
parameters, between central events and collisions with 
6 = 8 fm. In this case, K~ Y is then proportional to 
the system size R or, equivalently, to (1/5) dN/dy [see 
Eq. QJJ]. For small if -1 , we expect from Eq. Q that 
v%jt is proportional to K^ 1 (since v^ 5 " 110 scales like e), 
i.e., that v^/e scales like (1/5) dN/dy. Indeed, such a 
linear variation was found in Ref. Q (see Fig. 25), for 
both the RHIC Au-Au values and the SPS Pb-Pb values. 
It even turns out that both linear dependences match 
well, although the densities differ (the densities at SPS 
are typically 40% smaller than at RHIC), which suggests 
that the cross section (and c s ) used in estimating the 
Knudsen number can be taken as roughly constant across 
the different beam energies. 

We take the above-mentioned plot in Ref. 0] as evi- 
dence that local equilibrium is not achieved, since i^/e is 
steadily increasing with (1/5) dN/dy, without any hint 
at a saturation, for the whole accessible experimental 
range. A similar evidence comes from the variation of 
V2 with pseudorapidity 77 |l9l |20(: changing the latter 
induces a variation of the multiplicity density dN/drj, 
hence of the control parameter at fixed geometry, which 
is reflected in v 2 (q)- Finally, the decrease in the number 
of collisions as the transverse momentum increases also 
results in a departure from equilibrium seen on i>2 (pt) 
data and in transport calculations @- In the latter, 
for the bulk of the particles (jpr < 500 MeV/c), V2 already 
saturates for a of the order of 3 mb. For higher values 
of pt, the saturation occurs only above 10 mb. Surpris- 
ingly, however, the saturation value is significantly below 
that given by hydrodynamical computations, which calls 
for further investigation of the discrepancy. 

One can also study the effect of incomplete equilibra- 
tion on the fourth moment D4. Repeating the same ar- 
guments as for V2, one expects that the behaviour of v± 
as a function of K~ l is given, at least qualitatively, by 
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an equation similar to @. From the simple observa- 
tion that both u 2 an d V4 are proportional to K^ 1 for 
small K , one expects that v^jv\ decreases with K , 
reaching a minimum when the hydrodynamical regime is 
reached. Since it was shown in Ref. 21] that ideal fluid 
dynamics yields v±jv\ — 1/2 (a similar value was found 
in Ref. \ 22j within a specific hydrodynamical model), one 
expects the ratio to be larger than 1/2 if the system is 
not fully equilibrated, in agreement with the experimen- 
tal finding v±jv\ ~ 1.2 @|. 

Further predictions can be made to test the assump- 
tion of incomplete thcrmalization. In particular, study- 
ing smaller systems at the same center-of-mass energy, 
where the density is roughly the same (compare Tables [I] 
and[nj, one can obtain direct information on the depen- 
dence of i>2/e on the number of collisions per particle, 
K . Consider, e.g., Cu-Cu collisions at b = 5.5 fm, 
which corresponds to collisions with the same centrality 
as Au-Au collisions with b = 8 fm 16] . If the hydro- 
dynamical regime were reached, i>2/e would be roughly 
independent of the system size: from the values of the 
eccentricities in Table [HJ one concludes that the values 
of elliptic flow in Au-Au and Cu-Cu collisions would be 
then related by t>2(Cu) = 0.69 U2(Au). If, on the contrary, 
we are far from equilibrium, so that u 2 /e is proportional 
to K~ r , instead of constant, then the relationship reads 
Ua(Cu) = 0.34t>2(Au). This prediction can be completed 
by an analogous one concerning v^jv\: as argued be- 
fore, this quantity is a decreasing function of the number 
of collisions, hence should be larger in Cu-Cu collisions 
than it is in Au-Au collisions. Other predictions can be 
made, regarding the future experiments at LHC: if, as we 
have argued, the ideal-fluid limit is not reached even in 
the most central collisions at RHIC, then i^/e should fur- 
ther increase in Pb-Pb collisions at ^/snn = 5.5 TeV. In 
parallel, one can expect that v^jv\ will decrease, coming 
closer to 1/2. 

The case for early thermalization at RHIC rests on the 
argument that the vi measurements saturate the hydro- 
dynamical limit. However, the latter may well be under- 
estimated. Indeed, the pion data presented in Fig. 36 
of Ref. 2] suggest that the measured V2 in central colli- 



sions overshoots the value obtained by a hydrodynamical 
computation, V2(data) > V2 (hydro). We know from the 
above discussion that it is possible to increase the hydro- 
dynamical prediction for V2 by increasing the speed of 
sound. Now, in present hydrodynamical calculations, in- 
creasing c s means taking a harder equation of state, and 
this spoils the agreement with experimental momentum 
spectra, which require a soft equation of state |6J . It is im- 
portant to realize that this requirement is closely linked 
to the assumption that the system is fully thermalized 
(see |23j for an alternative to 3-dimensional thermaliza- 
tion) and, even more importantly, reaches chemical equi- 
librium: in a system in chemical equilibrium, there ex- 
ists a one-to-one relationship between energy per particle 
(probed by momentum spectra) and density (probed by 
dN/dy). Although fits to particle ratios by thermal mod- 
els are usually interpreted as evidence for chemical equi- 
librium, there is no direct experimental evidence that the 
particle density in the system obeys the laws of thermo- 
dynamics. In fact, as we have seen before, data indicate 
that kinetic equilibrium is not attained at RHIC, which in 
turn suggests that chemical equilibrium is not attained 
either. (Even recent hydro calculations, which assume 
early chemical freeze-out, start with a system in equi- 
librium |3.) If one drops the assumption of chemical 
equilibrium, energy per particle and density become inde- 
pendent variables and the transverse momentum spectra 
no longer constrain the equation of state. 

To summarize, relaxing the constraint of chemical equi- 
librium allows for a natural explanation of RHIC data on 
elliptic flow. Deviations from local equilibrium lead to a 
characteristic dependence of observables such as i> 2 /e and 
V4/V2 on the number of collisions, and this can be tested 
exp er iment ally. 



Acknowledgments 

We acknowledge the financial support from CEFIPRA, 
New Delhi, under its project no. 3104-3. 



[1] J.-Y. Ollitrault, Nucl. Phys. A 638, 195c (1998). 
[2] J. Adams et al. [STAR Collaboration], Phys. Rev. C 72, 
014904 (2005). 

[3] S. S. Adler et al. [PHENIX Collaboration], Phys. Rev. 

Lett. 91, 182301 (2003). 
[4] C. Alt et al. [NA49 Collaboration], Phys. Rev. C 68, 

034903 (2003). 
[5] J.-Y. Ollitrault, Phys. Rev. D 46 229 (1992). 
[6] P. F. Kolb and U. Heinz, in Quark Gluon Plasma III, 

World Scientific (Singapore), p. 634. 
[7] D. Molnar and M. Gyulassy, Nucl. Phys. A 697, 495 

(2002) [Erratum-ibid. A 703, 893 (2002)]. 
[8] D. Molnar and P. Huovinen, Phys. Rev. Lett. 94, 012302 

(2005). 



[9] D. Teaney, J. Lauret, and E. V. Sfmryak, 

nucl-th/01 10037 
[10] M. Gyulassy and L. McLerran, Nucl. Phys. A 750, 30 

(2005); E. V. Shuryak, Nucl. Phys. A 750, 64 (2005). 
[11] R. Baier, A. H. Mueller, D. Schiff, and D. T. Son, Phys. 

Lett. B 502, 51 (2001); ibid. 539, 46 (2002). 
[12] P. Arnold, J. Lenaghan, G. D. Moore, and L. G. Yaffe, 

Phys. Rev. Lett. 94, 072302 (2005). 
[13] A. Rebhan, P. Romatschke, and M. Strickland, Phys. 

Rev. Lett. 94, 102303 (2005). 
[14] D. Teaney, P hys. R ev. C 68, 034913 (2003) 
[15] P. Huovinen, [nucl-th/0505036 j 
[16] L. W. Chen and C. M. Ko, |nucl-tfi/0505044| 
[17] M. Knudsen, Ann. Phys. (Leipzig) 28, 75 (1909); E. M. 



5 



Lifshitz and L. P. Pitaevskii, Physical Kinetics, Perga- 

mon Press, (Oxford, 1981), page 48. 
[18] S. Gupta, hep-ph/0507210 
[19] T. Hirano, Phys. Rev. C 65, 011901 (2002). 
[20] U. W. Heinz and P. F. Kolb, J. Phys. G 30, S1229 (2004). 
[21] N. Borghini and J.- Y. Ollitrault, nucl-th/0506045 



[22] P. F. Kolb, Phys. Rev. C 68, 031902(R) (2003). 

[23] U. W. Heinz and S. M. H. Wong, Phys. Rev. C 66, 014907 

(2002) 

[24] T. Hirano and M. Gyulassy, nucl-th/0506049 



